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SHIMMING OF MRI SCANNER INVOLVING FAT SUPPRESSION AND/OR BLACK BLOOD 
PREPARATION 

The following relates to the diagnostic imaging arts. It finds particular application in 
magnetic resonance imaging, and will be described with particular reference thereto. However, it 
also finds more general application in nuclear magnetic resonance characterization techniques such 
as magnetic resonance spectroscopy. 
5 Magnetic field nonuniformities in magnetic resonance imaging scanners are known 

to introduce imaging and spectroscopy artifacts. These artifacts are particularly severe at higher 
magnetic fields and for certain techniques, such as for steady state free procession (SSFP) imaging 
in which magnetic field nonuniformities produce banding defects, and chemical shift imaging in 
which poor field uniformity produces spectral line broadening and shifting. 

10 Magnetic field nonuniformities are addressed in the first instance in design of the 

main magnet, which is typically a superconducting coil, although resistive main magnet coils are 
also used. Magnetic field uniformity is also improved through the use of selectively placed steel 
shims. Further improvement is obtained through the use of active shimming. In this technique, 
shim coils produce auxiliary magnetic field gradients that combine with the magnetic field 

15 produced by the main magnet to generate a more uniform overall field. Alternatively or 
additionally, offset currents can be applied to main gradient coils to provide active shimming. 
Advantageously, active shimming permits the radiologist to optimize field uniformity for a specific 
imaging session by performing preparatory pre-scan magnetic field measurements on the imaged 
region of the subject and adjusting the shim currents based on these subject measurements. 

20 Previous methods and pulse sequences for performing the magnetic field 

measurements and extracting shimming information therefrom have suffered from certain 
disadvantages. Accurate measurement of the magnetic fields has been hindered by artifacts 
introduced by resonance frequency shifts of fatty tissue, dephasing of the magnetic resonance in 
flowing blood, and cardiac cycling motion artifacts. Extraction of quantitative information has been 

25 further hindered by phase ambiguities across spatial regions of low magnetic resonance signal such 
as air in the lungs. As a consequence, active shimming has in practice employed collection of 
overdetermined data sets to ensure reasonable shimming accuracy. Even so, active shimming has 
generally in the past been limited to correction of first order magnetic field nonuniformities. 

The present invention contemplates an improved apparatus and method that 

30 overcomes the aforementioned limitations and others. 
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According to one aspect, a method is provided for shimming a magnetic resonance 
imaging scanner. At least one of fat suppression and moving blood suppression is performed to 
suppress magnetic resonance from at least one region of fat and proving blood First and second 

5 magnetic resonance echoes emanating from a generally columnar volume are measured. The 
measurements of the first and second magnetic resonance echoes are reconstructed into 
corresponding first and second generally columnar projection images. The first and second 
generally columnar projection images are combined to produce a magnetic fi^ld profile of the 
generally columnar volume. Selected magnetic field parameters are extracted from the magnetic 

10 field profile. Shim currents are computed from the selected magnetic field parameters. The shim 
currents are applied to magnetic field coils. 

According to another aspect, an apparatus is disclosed for shimming a magnetic 
resonance imaging scanner. A selective resonance suppression means is provided for performing at 
least one of suppression of magnetic resonance in fat and suppression of magnetic resonance in 

15 blood. A means is provided for measuring first and second magnetic resonance echoes emanating 
from a generally columnar volume. A means is provided for reconstructing the measurements of 
the first and second magnetic resonance echoes into corresponding first and second generally 
columnar projection images. A means is provided for combining the first and second generally 
columnar projection images to produce a magnetic field profile of the generally columnar volume. 

20 A means is provided for extracting selected magnetic field parameters from the magnetic field 
profile. A means is provided for computing a shim current from the selected magnetic field 
parameters. A means is provided for applying the shim current to a magnetic field gradient coil. 

One advantage resides in improved handling of phase ambiguities in phase-based 
magnetic field measurements. 

25 Another advantage resides in suppression of magnetic resonance from fat and 

blood regions during magnetic field mapping. 

Yet another advantage resides in more efficient active shimming by reducing the 
number of field projection measurements and avoiding overdetermined sets of linear equations in 
the shimming calculations. 

30 Numerous additional advantages and benefits will become apparent to those of 

ordinary skill in the art upon reading the following detailed description of the preferred 
embodiments. 

The invention may take form in various components and arrangements of 
components, and in various process operations and arrangements of process operations. The 
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drawings are only for the purpose of illustrating preferred embodiments and are not to be construed 
as limiting the invention. 

FIGURE 1 diagrammatically shows a magnetic resonance imaging system 
including active shimming. 
5 FIGURE 2 shows a preferred magnetic resonance sequence for performing the 

magnetic field measurements. 

FIGURE 3 diagrammatically shows details of the shimming processor of FIGURE 

1. 

10 

With reference to FIGURE 1, a magnetic resonance imaging scanner 10 includes main 
magnet coils 12, which are preferably superconducting coils, although resistive main magnet coils 
or permanent magnet material can also be employed. The main magnet coils 12 are energized to 
generate a substantially uniform main (B 0 ) magnetic field in an examination region 14. Magnetic 

15 . field gradient coils 16 produce gradients in selected spatial directions to spatially encode magnetic 
resonances that are generated by energizing a radio frequency coil 18. In FIGURE 1, a whole-body 
radio frequency coil 18 is shown; however, local coils such as head coils, phased radio frequency 
coil arrays, SENSE coils, and the like can be used instead of or in conjunction with the whole-body 
radio frequency coil 18 to excite magnetic resonances or to detect magnetic resonance echoes. 

20 A magnetic resonance sequence controller 30 coordinates and controls magnetic 

field gradient/shim controllers 32 coupled to the gradient coils 16 and a radio frequency transmitter 
34 coupled to the whole-body radio frequency coil 18 or another radio frequency coil to excite and 
spatially encode magnetic resonance echoes. One or more radio frequency receivers 36 coupled to 
the whole-body radio frequency coil 18 or another radio frequency coil detects, demodulates, and 

25 digitizes the magnetic resonance echoes and stores digital magnetic resonance samples in a k-space 
memory 40. A reconstruction processor 44 performs a Fourier transform-based image 
reconstruction or other type of image reconstruction to generate one or more reconstructed images 
from the stored k-space magnetic resonance samples. 

The reconstructed images are stored in an image memory 46, processed by a video 

30 processor 50 and displayed on a user interface 52, transmitted over a local computer network or the 
Internet, or otherwise processed. Preferably, the user interface 52 includes a display, printer, or 
other output device that allows a radiologist or other operator to view, render, or otherwise 
manipulate the reconstructed images. Moreover, the user interface 52 preferably enables the 
radiologist or other operator to communicate with the magnetic resonance sequence controller 30 to 
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create magnetic resonance imaging sequences, modify imaging sequences, execute imaging 
sequences, or otherwise control the magnetic resonance imaging scanner 10. 

Those skilled in the art recognize that the uniformity of the main (B 0 ) magnetic 
field is a substantial factor controlling the accuracy of the reconstructed images. B 0 field 

5 non-uniformities can introduce artifacts such as banding defects in steady state free procession 
(SSFP) imaging. To reduce such field non-uniformities, active shim coils 60 are selectively 
energized to produce auxiliary magnetic field gradients that compensate for non-uniformities of the 
main B 0 field. Additionally or alternatively, offset currents are applied to the main magnetic field 
gradient coils 16 to provide the active shimming. 

10 In order to determine appropriate shim or offset currents, non-uniformities of the 

main B 0 field are preferably measured using a preparatory pre-scan magnetic resonance sequence 
after the imaging subject is positioned in the scanner 10. For shimming of magnetic field 
non-uniformities of order N 62, a suitable preparatory magnetic resonance sequence is recalled via 
the user interface 52 from a shimming preparatory sequences memory 64. In a preferred 

15 embodiment, correction is performed at least for second order (N=2) magnetic field 
non-uniformities. The sequence controller 30 executes the recalled preparatory sequence to cause 
selected magnetic resonance echoes to be generated and measured, with the measurement samples 
being stored in a shimming k-space memory 70. A shimming processor 72 analyzes the measured 
echoes and computes suitable shim or offset currents to apply to the active shim coils 60 or the 

20 main magnetic field gradient coils 16 to compensate for the measured magnetic field 
non-uniformities. Preferably, the measurements are performed in the area to be imaged and with 
the imaging subject situated in the examination region 14. 

With reference to FIGURE 2, a preferred preparatory shimming magnetic 
resonance sequence 80 is described. The sequence 80 employs phase-sensitive imaging of a 

25 generally columnar volume defining a projection direction. Moreover, the shimming magnetic 
resonance sequence 80 provides suppression of magnetic resonance from fat tissues using a fat 
saturation preparation. This suppression is advantageous because the spectral resonance shift 
between fat and water causes apparent dephasing that interferes with phase-based shimming 
measurements. Similarly, the sequence 80 provides suppression of magnetic resonance from blood 

30 using a black blood preparation. This suppression is advantageous because rapid blood flow causes 
resonance dephasing that also interferes with phase measurements. 

In the preparatory magnetic resonance sequence 80, a black blood preparation 
includes a non-selective 180° radio frequency pulse 82 that inverts spins in the subject. A second 
180° radio frequency pulse 84 in conjunction with a slice-select gradient pulse 86 reinverts the 

35 static spins of the slice of interest. A delay following the two inversion pulses 82, 84 provides time 
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for blood to flow out of the slice of interest and for the magnetic resonance signal of the blood to 
relax toward the null condition. Before the blood is fully nulled, a third 1 80° radio frequency pulse 
90 in conjunction with a second slice-select gradient pulse 92 inverts the spins in the slice. A delay 
following the inversion pulse 90 is selected such that both inverted blood flowing into the slice and 
5 the magnetic resonance signal of the fat reach the null condition at a time when the slice is excited 
with a 90° radio frequency pulse 100 in conjunction with a slice-select gradient pulse 102. The 
magnetic resonance of the excited slice is refocused using a pair of 180° radio frequency pulses 
110, 112 in conjunction with two spatially transverse gradient pulses 114, 116 to define a generally 
columnar volume of interest 

10 Magnetic resonance echoes from the generally columnar volume of interest are 

measured with a short echo planar sequence 120 that reads a spin echo 122 and a field or gradient 
echo 124. The sequence 80 is a symmetric sequence: the refocusing 180° radio frequency pulses 
110, 112 are spaced at intervals T and 3T, respectively, from the 90° radio frequency excitation 
pulse 100, and the spin echo is read at about an interval of 4T from the 90° radio frequency 

15 excitation pulse 100. The field echo is read at an interval 4T+S from the 90° radio frequency 
excitation pulse 100, where the sub-interval S is selected to correspond to the multiple of the 
fat-water period. The fat-water period is the inverse of the frequency separation of the fat and water 
magnetic resonance signals. The sub-interval S is preferably a small multiple of the fat-water 
period, and more preferably the multiple is unity. With the spin echo and the field echo temporally 

20 separated by a multiple of the fat-water period, a phase difference between one-dimensional 
projection images reconstructed from the spin echo and the field echo, respectively, advantageously 
substantially corresponds to spin evolution under the influence of the local magnetic field. 
Contributions due to the fat-water chemical shift are advantageously suppressed by setting the 
sub-interval S to a small multiple of the fat-water period, as well as by the fat saturation 

25 preparation. 

The sequence 80 is preferably a symmetric sequence selected to provide 
compensation for slowly varying field terms due to longer duration eddy currents and field drift 
due to passive shims, booster iron, and the like. Phase shifts due to field terms with linear variation 
cancel at the second spin echo. Optionally, however, the preferred dual echo sequence 80 can be 
30 replaced by a single 90°-l 80° pulse pair. Moreover, while in the sequence 80 the fourth echo is read 
by the short echo planar sequence 120, those skilled in the art can readily modify the sequence 80 
to read flow-compensated third, fifth, sixth, or otherwise-indexed spin echoes. Still further, a 
stimulated echo can be substituted for the spin echo in the preparatory sequence 80. 

The sequence 80 is repeated for a plurality of spatial orientations of generally 
35 columnar volumes of interest. To extract magnetic field parameters of order N, P projections along 
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at least 2N+1 directions are measured. Thus, to extract second order magnetic field parameters 
where N=2, at least five orientations are measured For N^, at least seven orientations are 
measured. Measurement of a properly selected 2N+1 orientations of the generally columnar 
volume of interest allows for correction values for N(N+2) shim coils to be computed for magnetic 
5 field corrections up to and including order N. Optionally, more than 2N+1 orientations are 
measured; however, this produces an overdetermined system of linear equations and substantially 
reduces efficiency. Hence, in a preferred embodiment, the least number of orientations, namely a 
selected 2N+1 orientations, are measured. 

Using conventional spherical coordinate notation (6, <|>), appropriate orientations 
1 0 for the least number 2N+1 of generally columnar volumes of interest are selected according to: 

e * = ^TT^ 1+2k ) ' k = 0,l,...(N 9 -l) (1), 

and: 

15 

4=JL(l+2j) , j = 0,l,...(N,-l) (2), 

where Ne is a number of 9 coordinate values, N<|, is a number of § coordinate values. The 
product Ne'N^, is equal to the total number of collected projections P, and the value of N$ is chosen 

20 sothatN$>2N+l. 

As an example of the use of Equations (1) and (2), for second order shimming 
where N=2, the least number of orientations is 2N+1=5. Selecting Ne=l and N<|>=5, an 
appropriate set of orientations is given by Equations (1) and (2) as: (45°, 36°), (45°, 108°), 
(45°, 180°), (45°, 252°), and (45°, 324°). Equations (1) and (2) can be readily employed to compute 

25 appropriate sets of orientations for first order, third order, or higher order shimming as well. 

With reference to FIGURE 3, the shim k-space memory 70 stores k-space data for 
the spin echo and field echo acquired for each selected orientation of the generally columnar 
volume of interest. The shimming processor 72 processes this data to compute shimming currents 
or offset currents to implement the shimming. An iterative loop process 140 repeats the preparatory 

30 sequence 80 to acquire data for each orientation, and processes the k-space data for each orientation 
to compute magnetic field parameters along that orientation, which are stored in a magnetic field 
parameters memory 142. In a preferred embodiment, a bar reconstruction processor 150 applies a 
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filtered one-dimensional Fourier transform to reconstruct a spin echo columnar projection image 
152 from the spin echo data. The bar reconstruction processor 150 also applies the filtered Fourier 
transform to reconstruct a field echo columnar projection image 154 from the field echo data. 

A difference processor 160 combines the spin echo projection image 152 and the 
5 field echo projection image 154 to produce a shim bar projection 162 having a phase that is 
substantially due to spin evolution under the influence of the local magnetic field. Representing 
values of the two projections 152, 154 by complex numbers Z\ and Zz, and values of the shim bar 
projection 162 as complex number Z 3 , the shim bar projection 162 is suitably computed as 
Z 3 =Z2-Zi*, where Z x * is the complex conjugate of Zj. The shim bar projection 162 is a 

10 phase- wrapped shim bar projection. That is, the phase values of the shim bar projection 162 lie 
between 0 and 2%. A phase unwrap processor 166 unwraps the phase to remove phase 
discontinuities. Phase unwrapping is suitably accomplished by adding appropriate multiples of 2% 
to the phase values to remove phase discontinuities. 

Typical phase unwrapping algorithms can produce ambiguous or erroneous results 

15 in regions of low magnetic resonance signal. To address this difficulty, a region finder 170 
preferably analyzes the shim bar projection 162 to identify a largest usable region of substantial 
magnetic resonance signal. In one suitable approach, a global mean of the amplitude values of the 
phase-wrapped shim bar projection 162 is computed. Looping through the data, local mean values 
over short windows of typically three to seven data points are computed. The local mean values are 

20 compared with the global mean value. Transition points at which the local mean increases above 
the global mean indicate transitions between an unusable region and a usable region. Similarly, 
transition points at which the local mean decreases below the global mean indicate transitions 
between a usable region and an unusable region. The largest usable region is preferably unwrapped 
by the unwrapping processor 166, and this unwrapped shim bar portion is used in subsequent 

25 processing. 

The unwrapped shim bar projection or projection portion is fit to a polynomial of 
order N 62 by an Nth-order polynomial fitting processor 174 to extract the magnetic field 
parameters of order N 62 which are stored in the magnetic field parameters memory 142. In a 
preferred approach, the method of orthogonal polynomials is used in the fitting process; however, 
30 other polynomial fitting methods can be employed. The iterative loop process 140 invokes the bar 
reconstruction processor 150, the difference processor 160, the phase unwrapping processes 166, 
170, and the polynomial fitting processor 174 for each orientation selected using Equations (1) and 
(2), and the resulting magnetic field parameters for the various selected orientations are stored in 
the magnetic field parameters memory 142. 
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10 



A shim terms calculation processor 180 computes shim terms based on the 
magnetic field parameters of the various measured shim bar orientations. If shim orientations given 
by Equations (1) and (2) are used, then a least number of orientations are measured, and the linear 
equations relating the shim currents to the highest order magnetic field parameters are exactly 
determined. Advantageously, the choice of orientations according to Equations (1) and (2) ensures 
that the linear equations are fully decoupled and that the [W] matrix of Equation (4) is such that 
[W][W] T is a diagonal matrix even when an overdetermined measurement is made. In this case, 
computation of the shim currents reduces to algebraic calculations. The shim terms Ls Wm are given 



where W S him(0,<i>) denotes the angular part of a spherical harmonic function and A^r* denotes the 
fitted polynomial coefficients of the indicated order. On the other hand, if orientations different 
15 from those given by Equations (1) and (2) are measured it is possible to have an inconsistent set of 
linear equations relating magnetic field parameters and shim currents. Even if orientations are 
measured that lead to consistent sets of coupled linear equations, these may be overdetermined at 
the highest shim order or at least for the lower shim orders. In this case the shim terms calculation 
processor 180 can compute shim terms directly from the matrix equation: 

20 

[A] = [W][L] (4), 



where [A] is the polynomial coefficients matrix, [L] is the shim terms matrix, and [W] is the 

spherical harmonic functions matrix. Solution of matrix Equation (4) employs a less efficient 
25 algorithm that is suitable for solving overdetermined systems of linear equations. For 

overdetermined systems, singular value decomposition of the [W] matrix is preferably used to 

provide a solution that minimizes the effect of measurement error in the least squares sense. 

The invention has been described with reference to the preferred embodiments. 

Obviously, modifications and alterations will occur to others upon reading and understanding the 
30 preceding detailed description. It is intended that the invention be construed as including all such 

modifications and alterations insofar as they come within the scope of the appended claims or the 

equivalents thereof. 
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